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An artificial dielectric polymeric composite containing alloy particles which are oxide-coated for isolation and fer¬ 
romagnetic below room temperature for packing determination, is described. Induced magnetic dipole effects are 
calculated and compared with measurements of complex permeability in the centimeter wavelength range. Dielec¬ 
tric behavior, according to various theoretical model of heterogeneous media, is explored and correlated to experi¬ 
mental systems. Evidence of higher-order multipole effects and packing-limit effects in random composites is 
presented. Several advanced preparation techniques for fine powder were initialed.^. 
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SUMMARY 


The development of an interesting artificial dielectric is described. The metallic alloy particles 
(Ni 9 2 Cr g ) are oxide-coated for isolation and ferromagnetic below a Curie point of 165K to facilitate 
determination of volume loading in the composite. Particle size as well as volume loading are 
significant factors in the electromagnetic properties of the composites. Induced magnetic dipole (eddy 
current) effects depend markedly on the ratio of particle radius to skin depth which in turn depends on 
frequency and resistivity. The material and measurement parameters of this study span the range of in¬ 
teresting behavior. Measurements of the complex permeability agree satisfactorily with calculations. 
These phenomena have technological implications for real systems in extremes of frequency and/or 
temperature. An up-to-date exploration of theoretical descriptions of the dielectric properties of hetero¬ 
geneous media was made. They successfully account for the relatively simple Maxwell (Maxwell- 
Garnett. Clausius and Mossotti) behavior observed for artificial dielectrics of the present class. The 
effects of higher-order multipole interactions and geometrical packing limits calculated for regular lattice 
arrays are found in the random composites, with deviations that may reflect details of the composite 
microstructure. 

Work was also initiated on alternative passivation coalings and on the preparation of very fine 
powders by the spark erosion method. Such powders are intended for an empirical optimization of met¬ 
allic magneto-dielectric composite systems in the centimeter wavelength range. The spark erosion 
method presents several potential advantages of technological interest. 








I. INTRODUCTION 


This research is concerned with the development of a novel class of artificial dielectrics and a deeper 
understanding of its properties at millimeter wavelengths. Our starting point is the “Theoretical Radar 
Absorbing Materials Physics Model” developed during a classified program [1]. That study involved 
experiment and analysis of a broad data base obtained from measurements in the microwave range of 
frequencies (centimeter wavelengths). Many of the model features and predictions are, however, 
relevant to the millimeter wavelength region. Some of the more interesting implications of the work 
concern the dielectric property behavior of the composite system, one which may be regarded as a new 
class of oxide-passivated metal-insulator dielectrics. 

An overview of the research goals for the current three-year study may be stated in question form 
as follows: How valid is the RAM Physics Model in higher frequency regimes? Why does the dielec¬ 
tric behavior follow the simple Maxwell approximation, given the pitfalls that often beset metal- 
insulator composites? Has the artificial (magneto-)dielectric system been optimized in the centimeter 
wavelength regime? Can composites be fabricated which work well in the millimeter wavelength 
regime? How do they compare with other materials approaches? 

Alternatively, a more formal view of the program can be adopted from the Statement of Work given 
in the proposal. We reproduce as Table 1 the program schedule with time lines revised to fit the 
investigator's calendar, i.e., the contract start dale is 1 Aug 83 but formal signing and notification did 
not occur until about mid-September 1983. At an oral review/site visit by the ONR Scientific Officer 
on 12 June 1984, some task revisions were discussed and tentatively approved. These involved expan¬ 
sion of effort in Tasks 1 and 3, balanced by removal of Task 4. The rationale for these changes is 
developed herein. 

TABLE 1 

PROGRAM SCHEDULE 
AT START OF CONTRACT 





II. TASKS 1. 2, and 3 


The concurrent work on these three tasks during the reporting period has been so mutually interac¬ 
tive that it is better to consider their methods and results more or less together. The presentation of 
the segments will skip among the numbered Tasks in order to demonstrate a purposeful course. 

11.1 Development of an Artificial Dielectric 

With initial emphasis of the program on the dielectric aspects, an important phase has been the 
development of a suitable laboratory artificial dielectric. The model idea of a dielectric which consists 
of a series of metallic globules separated from each other by insulating material goes back at least to 
Faraday in 1837. An excellent history of inhomogeneous media concepts has been recorded by Lan- 
dauer [2], 

11.1.1 Criteria 

Our criteria are that the metal (alloy) powder be non-magnetic at room temperature but ferromag¬ 
netic at low temperature to facilitate determination of the volume loading, p. The latter quantity is key 
to quantifying the dielectric behavior. We also require that the powder be capable of developing a thin 
insulating oxide layer on each particle. We contrast this with some other work in which researchers 
have encouraged formation of a thin oxide coating in powder preparation to discourage initial 
agglomeration but kept it thin enough to allow metal-io-metal contact between the particles under 
sufficient pressure [3]. 

The choice of candidate alloys thus involved magnetic and metallurgical consideration. Added to 
these is the challenge of fine metal powder preparation. For this task area we made use of our gas- 
atomization facilities. In this process [4] molten alloy is forced through a specially designed nozzle and 
the liquid stream is disrupted by high velocity gas blasts to create fine droplets. These solidify as they 
descend in a large volume settling chamber. In one of our facilities liquid substrate quenching is a spe¬ 
cial feature, i.e., a blast of water spray follows the gas jets. 

11.1.2 Consideration of Candidate Metallic Systems 

Desirable alloy systems were therefore such that one of the elemental components would tend to 
form a stable and at least partially protective oxide. Some usual examples of this category are Al. Cr 
and Y. Folding these in with the magnetic considerations, yielded several suggestions, e.g.. Ni with Cr. 
Gd with Y, and (Cu-Ni) with Al. The first two systems are reasonably well known magnetically in 
handbooks [5] or in more recent (6] research while the third has had little exploration [7], The Gd-Y 
combination was quickly rejected (by our metallurgical colleagues) as explosively undesirable. 
Although the solubility of Al in the Ni-Cu host is modest, at about 10 at % (—5 wt %), there appeared 
to be a fair prospect of a Curie point below room temperature. We chose to prepare an exploratory 
ingot with the composition (by weight) Ni90Cu5AI5. Finally for Ni-Cr, the wide range of solubility and 
well-known magnetic properties enabled us to specify “with confidence" an ingot composition of 
Ni92Cr8 with every expection of a Curie point near 240K, and a low temperature saturation magnetiza¬ 
tion of about fr oo (0°K) — 21.8 emu/g. 

11.1.3 Magnetic Characterization of Candidate Metallic Systems 

Thermomagnetic scans, magnetization vs. temperature at 50 Oe, along with low temperature (6K. 
20K) saturation magnetization data, were measured. For this purpose we have a vibrating sample mag¬ 
netometer system with capabilities from 4K to 1370K (—10000 and fields up to 25kOe. The scans are 
shown in Figure 1. As for the Ni-Cu-AI ternary, its T c centered at 256K but there was a risk of some 
residual undesirable ferromagnetism at room temperature. The opposite problem presented itself with 
Ni92Cr8. Both the Curie point and the magnetization fell short of confident expectations by a factor of 
2! In searching for the ’'fault” here, we verified the chemical composition of the ingot, and also 
prepared a second ingot with composition Ni94Cr6. Its T c , measured as 237K. again compares poorly 
with the literature value of 335K! At this point we were facing an unexpected, annoying and delaying 
problem. A resolution of the dilemma was clearly warranted because this Ni-Cr system remained the 
most attractive of the ones considered as candidates. 
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Figure 1. Thermomagnetic analysis of candidate alloys and atomized powder. 


II. 1.3.1 Exploration of Magnetic Characterization Problems 

The solution of the problem and the “validation” of our own work came from a detailed look at the 
published literature both backward from the handbook data and forward to the present time via a library 
bibliographic search. The magnetic measurements cited in the standard References [5] are based on 
experimental work from the Strasbourg laboratory directed by P Weiss, a world-famous organization 
that set most of the standards for current magnetics research. In studying the primary sources, the 
research of Marian [8] and Sadron [9], it emerges that they used identical alloy samples which had been 
prepared earlier by Safranek [10], also at Strasbourg. In the preparation, he had been concerned with 
the purity of the available starting materials, especially the chromium. Indeed, he ultimately prepared 
his own Cr by an electrolytic process, to combine with high purity Mond nickel. Satisfactory metallur¬ 
gical micrographs were shown, indicating single phase alloys. Inasmuch as our Curie points are 
markedly lower, a suspicious investigator may still want to raise the purity question. How much resid¬ 
ual Fe or Co in either component would be needed to shift the results into “accord"'.’ More tolerance 
is in order with regard to the magnetic saturation because the French workers had to extrapolate down¬ 
ward from 110K toward absolute zero. However, their results are 40% to 100"'!. greater than ours! 


As one looks forward from the standard references into the most recent two decades, a different pic¬ 
ture emerges. A group at Leiden [11] in 1962 noticed some differences in saturation magnetization for 
alloys of 5.6 at % and 9.1 at %. and in an estimate of T c for the former alloy, but politely called these 












"small discrepancies." Next, a rather thorough investigation in 1972, at Strasbourg by Besnus et 
al., [12] (on newly prepared samples and including low-temperatures) presented results which differed 
quite a bit from the early data though not quite as much as the Leiden results. They diplomatically 
omit to mention that they differ from their colleagues of an earlier generation! The emerging picture is 
displayed in Figure 2 for Curie points and in Figure 3 for saturation magnetization. (Compositions are 
in atom percent, as used by these research workers. In the region of interest to us, i.e., 6-8 wt%, the 
amount in atom % is about one integer higher). More recently a study from Japan [13] gave results in 
close agreement with Besnus et al., and did make a significant graphical display of the differences with 
the early work. Lastly, Simpson and Smith in 1982 report [14] especially on compositions with quite 
low Curie points, in close agreement with our measurements. 

While examination of the comparative results in Figures 2 and 3 does not fully settle the question of 
which set is correct, the existence of a major discrepancy with the early work is established. The 
present work on the several ingots sets the base line for subsequent magnetic evaluation of powder 
properties and ultimately, for our original goal, of volume loadings in artificial dielectric composites. 



Figure 2. Ferromagnetic Curie points as a function of composition in Ni-Cr solid solution 
alloys including work of other investigators as well as present results. 







Figure 3. Saturation magnetization per average atom vs. composition for Ni-Cr alloys in¬ 
cluding available literature data. 


11.1.4 Preparation and Magnetic Characterization of Metal Powder 

At this stage we were ready for powder preparation. We anticipated (correctly) that there would he 
some composition shift in this step toward a higher Curie point and thus we selected the ingot Ni92C’r8 
(about 9 at. % Cr) in order to keep the alloy non-magnetic at room temperature. 

Two atomization trials of the Ni92-Cr8 alloy were performed. In both cases the melt was atomized 
while at a temperature of 1870K by argon gas at 4.28 MPa. The high melt temperature was employed 
to reduce the melt surface tension and melt viscosity while the high gas pressure resulted in a high gas 
jet velocity. (Mach Number =3.2). Previous analysis has shown that fine powder sizes are favored by 
increased gas velocity and decreased melt surface tension and melt viscosity 

The two trials were performed using two different versions of the GE-CRD-developed close-coupled 
gas nozzle. The second run uses an improved form of the close-coupled concept, designed to provide 
increased yields of fine powder. 

The results of the two runs are shown in Figure 4. The particle size distributions between 3 ~n ni 
and 250/am were determined by sieving. The particle size distributions below 37 M m were obtained b> 
extrapolation assuming the curves to be log-normal. (This assumption has been verified in other 
work). 








The thermomagnetic scan of this powder (sieved to <37/um) is shown in Figure 1. An average 
Curie point of 165K results, some 68K upward from the ingot. It is broader, suggesting a range of T c 
values. 

Referring to Figure 2. this T c translates in composition to about 8 at. %Cr or 7.15± 0.2 wt"/hCr. Us¬ 
ing this composition, there is satisfactory accord of the saturation magnetization “line" in Figure 3 with 
the "present work" data. 

For the practical record, the magnetization values at low temperature (6K or 20K), extrapolated to 
infinite field, are 21.2 emu/g and 12.0 emu/g for the ingots Ni94Cr6 and Ni92Cr8 respectively. For the 
powder we find 16.2 emu/g. 

II.1.5 Evaluation of Surface Condition of Metal Powder 

The shift in composition during powder preparation corresponds to a net depletion of Cr in the alloy 
particle. It is tempting to attribute this to an oxide coating of Cr 2 0 3 . It is much more difficult to exam¬ 
ine the surface coating of spherical particles by microprobe techniques, than the surface scales on bulk 
flat samples. We are, however, attempting to do this on these powder spheres. Conventional oxidation 
studies on the Ni-Cr system, albeit under isothermal conditions, indicate a moderately complex scaling 
with NiO, Cr 2 0 3 , Ni-Cr spinel oxide and other components [151. 
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If we naively (and tentatively) hold to a view that the coating may be Cr 2 0 3 , we can estimate its 
thickness from a mass balance. Using a mean diameter of about 23/xm we calculate a shell coating 
whose thickness is about 0.1/un, which is not unreasonable. A sketch of this rough calculation is as 
follows: 

[(w/o Cr in ingot) x vol. sphere x p alloy ] = 

= l(w/o Cr in powder) x vol. sphere x p alloy ] + 

+ [(w/o Cr in oxide) x (vol. oxide shell) x p oxide ), 
vol. sphere = 47rr 3 /3 ; vol. shell = 4irr 2 Ar, 

(w/o Cr in ingot - w/o Cr in powder) x <diam.> x p 3) | 0y 
(w/o Cr in oxide) x 3 x 2 x p 0Xlde 

We have some results from microprobe techniques on the Ni-Cr powder. As-prepared powder was 
pressed onto and into an aluminum substrate block for examination by Auger electron spectroscopy. 
The count rates for the various elements, Ni, Cr and O, were monitored as a function of sputter time. 
The latter is proportional to depth into the particle being examined, as material is sputtered away. The 
probe spot used was smaller than the particle diameter (from the powder size fraction +20, - 37pm). 
Count rates are expressed in atomic percent, after appropriate scaling in the instrument. [Physical Elec¬ 
tronics Model 600 Auger Microprobe]. The “depth” probing may not be completely simple because 
the pressure used in preparing the specimen may have deformed the spherical alloy particle. As an 
alternative, we hope to obtain a profile on a polished particle as the probe traverses the edge coating. 

The depth profile by Auger spectroscopy is shown in Figure 5. While not the simple naive picture 
we hypothesized, it does have features one might expect. At the surface (t = 0, sputter time) there is 
a major build-up of oxygen , from a residual low level representing a sort of background; and a modest, 
~50%, build-up of chromium from about 12 at.% to about 18 at .%. Concurrently, there is a major 
decrease of nickel, albeit not to the near-zero-level hypothesized. We can attempt to rationalize these 
data in terms of NiO and Cr 2 0 3 , there being no evidence for Ni metal in the magnetic data. (Oxidation 
studies of this system [15] always assume that Cr goes to Cr 2 0 3 , or more complex spinels or solid solu¬ 
tion). The "surface composition,” Ni 032 Cr 0 .ig O 0}0 , is roughly equivalent to 3Ni01 Cr 2 0 3 , which 
corresponds to Ni 0 2 7 Cr 0 ,i 8 0 o54 . 

If we oversimplify the surface as a uniform shell (again), we can invoke a mass or atom balance to 
determine how much material is in the shell. Using the chemically and magnetically determined alloy 
compositions, we write 

Nio.91 Cr 0 .o9 + 6zO— 1 • y Ni 092 Cr 0 0 8 + z (3Ni01 Cr 2 0 3 ). 

Balance equations for Ni and Cr yield two simultaneous equations in y and z, which are readily solved 
to yield y = 0.969, z = 0.00625. Approximately 2% of the original nickel goes into the oxide shell, 
while about 14% of the original chromium makes the transition. This does confirm the general stronger 
tendency for Cr to oxidize. Expressing the chemical balance in terms of the mass of one average 
(23pm) particle and the number of moles of alloy in that particle, we can obtain the masses of NiO and 
of Cr 2 0 3 . Converting to volumes (via their separate specific gravity values), and equating their total 
volume to a "shell” volume, we arrive at the oversimplified thickness of shell of 0.22pm. This is a 
factor of two greater than the more naive estimate, but it too is only a rough (albeit reasonable) ap¬ 
proximation. 

II. 1.5.1 Resistance Tests of Metal Powder Compacts 

As a final practical acceptance test for the alloy powder, we attempted to verify that the powder par¬ 
ticles were electrically insulated from each other. Conceptually, a check of resistance in powder com¬ 
pacts appears simple. In practice, it is less so. A simple device was constructed consisting of solid 
copper electrodes which can slide into a plastic cylinder, 1.25 cm inner diameter. We used this with a 
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SPUTTER TIME (min.) 


Figure 5. “Depth” profile by Auger electron spectroscopy of Ni-Cr powder compact: Com¬ 
position vs. sputter time (proportional to depth). 

digital multimeter capable of measuring up to 20Mft resistance. A small amount of alloy powder in 
the cylinder gave erratic results, dependent on pressure, time, powder column height, etc. 
Confirmation of this behavior was found in the literature. Ottavi, et al„ [16] describe a conceptually 
elegant experiment to study percolation concepts, i.e., at what composition do continuous electrical 
paths develop when various mixtures of conducting and insulating spheres are studied? Their experi¬ 
ments were carried out with monodisperse 5mm diameter plastic spheres, a fraction of which had been 
electroplated with a thin copper coating. It is helpful to quote from their paper, as follows: “In order 
to obtain stable conductance measurements we find it necessary to apply pressure to the piston bearing 
the upper electrode. This lack of stability is not that due to fluctuations introduced by the randomness 
of the mixture. Indeed, even if the system is reduced to only two spheres, the measured resistance is 
not constant in time, or reproducible between trials, unless a sufficient force (about 10 kgf) is applied. 
— We observe, moreover, a hysteresis phenomenon 

Following these hints, we obtained roughly reproducible behavior. Powder column height was kept 
to 1.0 to 1.2 mm. Significant arm/body pressure, with some electrode piston rotation, was employed. 
In addition, a calibration scheme with other variously prepared powders was introduced. These 
included a hydrogen-reduced ferromagnetic metal powder, and a deliberately well-coated magnetic met¬ 
allic powder. Under these conditions, we obtained the results given in Table 2. On the basis of these 
results, we conclude that the gas-atomized Ni-Cr powder prepared from the ingot Ni92Cr8 should meet 
the insulated particle criterion for our artificial dielectric. 
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Table 2 

RESISTANCE OF ALLOY POWDER COMPACTS 


Sam pie-Powder 

Approx. Resistance 

Ni-Cr < 37/am 

Reduced Magnetic Metal < 37/am 

Well-Coated Magnetic Metal 14 < d < 27/am 

i o 5 — i o 6 n 

io 1 n 

> io 7 n 


II.1.6 Preparation of Artificial Dielectric Composites 
U. 1.6.1 Selection of Binder for Artificial Dieletrics 

The goals of the program call for investigation of artificial dielectrics in two different insulating 
binders. Our plans have been to use one polymeric binder (for which «' is typically 3±0.5) and one 
inorganic binder for which <' may be in the range 5 to 10. For the first of these there is an easy choice 
for a laboratory system, i.e., polyurethane, known commercially as Permuthane™ (Division of Beatrice 
Foods, Corwin St., Peabody, MA, 01960). It has been used extensively in composites and our colla¬ 
borators at GE-Aircraft Engine Group, Evendale, Ohio have carried out the preparation of test speci¬ 
mens for us. For the second binder, we asked them to try a proprietary developmental inorganic 
binder. This has not worked out well, as will be described below. We are currently discussing the pos¬ 
sibility of using another moderately high dielectric constant proprietary binder as a substitute. 

H. 1.6.2 Preparation of Polyurethane Test Speciments 

For our initial series of measurements, a set of nine test coupons was prepared. These involved 
several particle sizes, for reasons which will be explained in a later section. The basic specifications are 
given in Table 3, including measured values of p determined as described below in Section 11.1.6.3. 

The preparation followed guidelines for curing the polymer as prescribed by the manufacturer. Suit¬ 
ably chosen amounts of alloy powder and binder powder were combined and blended for several hours. 
The mixture was leveled into a 2.5 cm x 5 cm (1 inch x 2 inch) die to a thickness such that the final 
coupon would be ~1 mm (—0.040 in.). The loaded cold die was placed on hot press platens, heated 
under light pressure from both platens, then with increased pressure, reaching a curing temperature of 
177C where it was held for 5 min. It was subsequently cooled under cold running water. 

The weighing calculations were based on a cured Permuthane density of 1.06 to 1.11 g/cm 3 . (For 
sample ONR 1, we determined 1.11 g/cm 3 .) These densities, as well as those for the composites are 
obtained from dimensions on regular shapes and weight. This avoids the possible problems of porous 
media, etc., but does permit small uncertainties in the determinations. For the alloy, we used a density 
of S.79 g/cm 3 based on the magnetic determination of composition and an x-ray lattice parameter value 
of a 0 = 3.529± 0.006A. 

11.1.6.3 Evaluation of Volume Loading 

As indicated in Section II.1.1 we place emphasis on a straightforward and non-destructive approach 
to determining the volume loading, p, of metal in the artificial dielectrics, as a guide to analysis of their 
electromagnetic properties. Our approach to this is through magnetic measurements at low tempera¬ 
ture. having chosen an alloy metal that has its Curie point below room temperature. 

A few simple formulas are helpful here. The subscripts a, b. c and h will refer to alloy, binder, 
composite and holes or voids respectively. Although we use the symbol p to designate volume loading 
in the bulk of this report, for this discussion we shall use the symbol f with the superscript w or v to 
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Table 3 


DESCRIPTION OF POLYURETHANE 
TEST COUPONS (ALLOY-Ni, 2 Cr 8 ) 


Sample # 

Alloy Particle 
Size 
(in fim) 

Nominal Alloy 
Volume 
Loading, p 

Measured 
Volume 
Loading, p 

ONR 1 

- 

0 

- 

ONR 2 

+20 - 37 

0.40 

0.413 

ONR 3 

+ 20 - 37 

0.30 

0.307 

ONR 4 

+ 20 - 37 

0.20 

0.200 

ONR 5 

+ 20 - 37 

0.10 

0.112 

ONR 6 

+ 10-20 

0.40 

0.396 

ONR 7 

-10 

0.40 

6.384 

ONR 8 

-10 

0.20 

0.208 


indicate weight fraction or volume fraction. The symbols <r, M and p designate magnetization per 
gram, magnetization per unit volume (e.g., cm 3 ) and specific gravity (or density in g/cm 3 ). 

For specific magnetization, or magnetization per gram, of composite we write 
<r c - zr a f a w + cr b f b w + tr„ f h w , 

but obviously cr h - 0 and f" - 0 at all times, and thus f a w + f b « 1. Therefore we have <r c * <r a f* + 
cr b (l-f a w ). This expression become trivial but useful when one notes that cr b = 0 in most cases of 
interest. Therefore the formula for weight fraction of alloy becomes f a * - cr c /<r a . Simply the measure¬ 
ment of or c along with a knowledge (or independent measurement) of the alloy magnetization gives its 
weight fraction. 

For volume magnetization and the related volume fractions, we obtain in the same manner 
M c - M a f a v + M b f b v + M h f h v , 

but this time f b may be different from zero although M h = 0 at all times and M b - 0 in most cases of 
interest. Note also that f a + f b + fh ” 1, by definition. The formula for volume fraction of alloy, p, is 
thus 

p = f v = Me = Pc q~c =fW Pc 
a M a p a cr a a p a 

This requires the additional measurement (or knowledge) of the density of the composite and of the 
alloy. With these measured values, on whole samples, the results are on firm ground and contain no 
approximations or idealizations. We discuss determination of the void fraction in a later section. 

As for the case of the binder alone, we have obtained more consistent results in determining density 
in the “primitive” way from dimensions on regular shapes and weight. For many electromagnetic con¬ 
siderations it is the overall volume property which is of interest, so that this simple method which 
includes voids, cracks, etc., is the proper one. The results for the volume loading determinations are 
given in Table 3. The agreement between nominal and measured values is very satisfactory and should 
be characteristic of polymeric-type binders when appropriate care is taken. With other categories of 
binders such agreement cannot be taken for granted. 
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II. 1.6.4 Preparation of Inorganic Binder Composites 

For a second binder choice for the composite artificial dielectrics, we tried to use a proprietary inor¬ 
ganic binder currently under development at GE-Evendale. It is intended for a similarly-sized but 
chemically different dispersoid (loading material). Typically, with inorganic binders there is a high tem¬ 
perature cycle to bring the binder to its final state. In our past experience [1] such a high temperature 
excursion may be accompanied by chemical reaction between binder and dispersoid particle. In the 
present case we requested 4 composites, p = 0.1, 0.2, 0.3 and 0.4, each with the largest particle size 
+ 20-37/xm of the Ni$ 2 Cr 8 powder. In each case, the resulting composite was ferromagnetic at room 
temperature. Indeed a thermomagnetic scan showed Curie points up to 350C (pure Ni, T c = 3580 
with strong suggestions of composition variations within the magnetic phase: 

Under these conditions it is very difficult to know the alloy composition in the composite. There¬ 
fore we could not determine volume loading without destruction and visual or probe analysis of the 
composite. 

Still another proprietary binder with moderately higher e' is being developed at GE-Evendale. (e = 
5). In the coming period we hope to be able to make composites with this binder which is less 
demanding of the dispersoid. 

Therefore, the sought-for comparison with different permittivity binders has been delayed. We do 
have a few data points from other concurrent work with this binder which support the expected trend. 

II.2 Magnetic Dipole and Particle Size Effects 

At first consideration it might seem that purely electric dipole phenomena would fully describe the 
high frequency response of an artificial dielectric composed of strictly non-magnetic materials (in con¬ 
trast to prior work on artificial magneto-dielectrics ID). Indeed this appears to have been the initial 
point of departure for analysis over the past rather active decade of research on selective absorbers of 
solar and infrared radiation [17,18}. However, it was quickly realized lhat, if the diameter of the con¬ 
ducting particle is not negligible compared to the skin depth, there is a second, rather classical, mecha¬ 
nism for absorption arising from the effects of the magnetic field of the incident radiation. In the time- 
varying magnetic field, there will be an induced eddy current magnetic polarization and an eddy current 
loss. The effective permeability of the composite artificial dielectric will thus not be unity even though 
the d.c. permeability of the constituents may be exactly unity. These magnetic effects are size- 
dependent while the electric dipole polarization and loss phenomena are not. Indeed the magnetic con¬ 
tribution dominates for metallic particles above about 50A in size (17,19)! 

While this feature was not spelled out in our initial plans for this research program, calculations 
presented below show that much interesting behavior occurs in relevant frequency and particle size 
ranges. Therefore, we have put emphasis on obtaining and measuring different powder size fractions. 
This is the expansion of scope referred to in Section I above, which we have requested to be balanced 
by suppression of Task 4 of the original program schedule (Table 1). The change has also caused us to 
subcontract out a small effort at sizing some of the alloy powder by air-classification. Sieving is a sensi¬ 
ble mode only for small quantities. This subcontract was not envisaged when the original program 
plans were submitted. 

II.2.1 Magnetic Polarizability of a Single Particle 

The magnetic response of an isolated isotropic conducting sphere has been discussed nicely by Lan¬ 
dau and Lifshitz (20). Following their notation. 

M - a H 

where M is the magnetization (magnetic moment per unit volume), H is the external magnetic field 
and a is the (complex) magnetic polarizability. For a sphere they derive (with or — a' + i a"), 

3 8 [sinh (2a/8) - sin (2a/8)] 

2 a (cosh (2a/8) - cos (2a/8)) 
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Figure 6. Magnetic polarizability of a conducting sphere. Negative real part, -a', vs. the 
dimensionless ratio of sphere radius to skin depth. 


MAGNETIC POLARIZABILITY (REAL PART) 
CONDUCTING SPHERE (-«') 







In Figure 7, the imaginary part of the magnetic polarizability is presented, again as a function of a/8. 
This quantity is intrinsically positive. It vanishes in both of the extremes of small or large radius (or 
frequency); in the former the limiting behavior is as (a/8) 2 , while in the latter it decreases in propor¬ 
tion to 8/a. The maximum loss occurs at a/8 - 2.3. while remaining above half its maximum for the 
(a/8) range of 1.2 to 7. 

II.2.2 From Single Conducting Sphere to Composite 

The magnetic polarizability formulas given in the preceding section pertain only to the spherical 
shape and the externally applied magnetic field. In describing magnetic properties in general, as well as 
in the present case of a dispersion of small magnetic particles, one must proceed carefully in accounting 
for the various demagnetizing effects, internal and external fields, etc. For example, in the present case 
one may have to distinguish between the demagnetizing coefficients of the overall composite sample, of 
the individual particle, and of the cavity “constructed” around the particle in analyzing for interactions 
between the particle’s dipole moment and the magnetization of the particles. In our model, the latter 
interactions are considered in the mean-field approximation which is accurate in the dilute limit. (We 
shall describe some of its weaknesses in a later section when we deal with the dielectric behavior). As 
another example, particularly relevant to conventional magnetic materials, i.e., ferromagnetic ones, the 
quantity of fundamental importance is the intrinsic or internal property, be it susceptibility, permeability 
or polarizability. 

For the artificial dielectric composites studied here, we have the simple case of spherical particles 
(demagnetizing coefficient = 4fl73), spherical cavity (same) and overall sample configurations for 
which we can ignore demagnetization. The latter come about because the measurements are carried out 
with closed flux configurations (toroid in a co-axial line) or thin slabs with the incident H-field parallel 
to the slab plane. 



a/8, RADIUS/SKIN DEPTH 

Figure 7. Magnetic polarizability (imaginary or loss part) of conducting sphere in time- 
varying field, a" vs. (a/8). 
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Having hoisted the flags of caution, we shall not repeat the derivations which may be found in 
appropriate textbooks, but shall merely present some key formulas. The intrinsic (or internal) mag¬ 
netic susceptibility of the sphere (relative to the internal field) is related to the polarizability derived 
above as: 


For a composite of spherical particles with volume loading, p, the susceptibility relative to an externally 
applied field is: 


where N s is the demagnetizing coefficient for the overall composite sample shape. In the common case 
herein, we have N s = 0, and thus we obtain 


l + ^d-p) Xi 

Upon substituting the relation between and a, this formula becomes 


This appears so simple that one presumes it could have been found more directly. However, we feel 
the excursion into details and warnings is appropriate. Finally for the permeability, we have, by 
definition, Mexi * l+ 4 ?r Xexi and then 4 


At this point we take note that a and hence n are complex, with the forms a « 
Making these substitutions and continuing the manipulation, we obtain 


-jj l + yp«'- 2 (y— p ) 2 [(a ') 2 + (a") 2 l , 


fn Figures 8 and 9 we show the behavior of ji ex , and for values of p from 0.1 to 0.5. As will be 
discussed later, the empirical limit for random packing of monodisperse spheres is about 0.63, but in 
practice one rarely goes much beyond about 0 5. In Figure 8 , we note that the negative real polarizabil¬ 
ity manifests itself in a substantial lowering (diamagnetic tendency) of n' from unity. In Figure 9. the 
family of curves for fi" are straightforward reflections of the behavior of a” (as given in the formula) 
These formulas and graphs are in a universal dimensionless form, expressed in terms of a/ 8 . It 
remains to demonstrate that the changing properties are relevant to the materials and electromagnetic 
parameters of this investigation. 

II.2.3 From Model Composite to Practical Application 

At this point we transform the model calculations, expressed in dimensionless ratios of material 
parameters, into real units that are measurable. The major parameter is the resistivity, followed by par¬ 
ticle sizes and frequencies of interest. The resistivity. p EL , for the alloy Ni 92 Cr* is 55 x 10~ 8 llm [211. 
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Figure 8. Magnetic permeability (real part) calculated for composites of conducting spheres 
at various volume loadings, p; ft’„, vs. (a/8). 


When this is substituted into the formula for skin depth, 8. we obtain at 10GHz, 8 = 3.7/um. Proceed¬ 
ing a step further with the consideration that peak a" (and n”) occurs at a/8 = 2.3, we find that 
•'important” sphere diameters , d c are as follows: 

d c = 17 jun at 10GHz 

d c = 7.6fitr\ at 50 GHz 
d c = 5An at 100 GHz 

Obviously there is relevance for a powder with particles below 37^m (-400 mesh). 

We show in Figures 10 and 11, the calculated model behavior for n' and n" for this alloy composi¬ 
tion, for a volume loading, p = 0.4, and for various sphere diameters of interest, i.e., 5, 10, 20 and 
30/*m (Similar calculations for other p-values are presented below). It should be noted that our calcu¬ 
lations above have been in the c.g.s. system for which n', m" are dimensionless. As we move closer to 
real electrical measurements in the “more practical” m.k.s. SI system, these expressions or graphs of 
n n’ (and later of «',«”) are to be considered as relative to the permeability (or permittivity) of free 
space. In other words we have nr ” m’/mo. etc. 

The qualitative remarks made for Figures 8 and 9 remain true for these figures. However, the 
separation of the features for the various sphere diameters is rather dramatic. It certainly invites exper¬ 
imental demonstration in that the range of sizes and of frequencies for this investigation spans the cal¬ 
culated behavior very nicely. It is for these reasons that we prepared samples with different particle 
sizes (and volume loadings) as listed in Table 3; i.e., -10/am, + 10-20/um, and + 20-37//m. 

It is instructive to extend the critical diameter/frequency consideration somewhat further on either 
end. If we move to the far-infrared, to a frequency of 1 THz (1 x IO ,2 Hz) alternatively expressed in 
wave numbers as 33cm" 1 , the critical diameter for maximum absorption in this alloy would be 1.7/xm 
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(a/8) RADIUS/SKIN DEPTH 

Figure 9. Magnetic permeability (imaginary part) calculated for composites of conducting 
spheres; vs. (a/5) for various volume loadings, p. 

which is in line with predictions for such research [19]. This domain is of interest for selective 
absorbers of solar or i.r. radiation. At the other extreme, that of low frequency, the critical diameter at 
10 kHz is 1.7 cm. This is in the regime for induction heating! 

In each of these two limiting regions one is unlikely to find results obtained and presented in terms 
of the complex permeability (and permittivity). Instead the focus is on absorption per se, because of 
measurement techniques as well as because of practical interest. For these reasons, measurements of 
the permeability reported here make a useful contribution to the understanding of artificial dielectric 
materials in the microwave frequency region and they relate to the nearby far-infrared region also. At 
present, we are aware of only one report which presents these kinds of measurements and analysis. It 
deals with systems of millimeter-sized ferromagnetic steel spheres over the kilohertz and megahertz fre¬ 
quency range [22], 

There is yet another interest in such model results, a technological one. The RAM physics 
model [1] predicts that high frequency loss in a magneto-dielectric should become negligible beyond 
certain frequencies and/or temperatures (depending on material parameters). When the induced mag¬ 
netic effects are included, one may expect to find lossy properties outside the limits indicated previous¬ 
ly. An example of this is presented in the classified Appendix. 











Figure 10. Relative real magnetic permeability, /* r ', vs. frequency calculated for composites 
of conducting spheres of various diameters; composite volume loading, p = 0.4; 
sphere alloy resistivity, p E l - 55 x 10 ® II m. 


II.2.4 Measurement of High Frequency Magnetic Permeability and Permittivity 

Measurements of the constitutive parameters, complex p r and complex « r . are being made in the 
centimeter wavelength range, —3-18 GHz. at GE-CRD. The equipment is based on a Hewlett-Packard 
Network Analyzer. The major building blocks of the system are the 8620A Sweep Generator, the 
8410A Network Analyzer, the 8411A Harmonic Frequency Converter, the 8412A Phase-Magnitude 
Display, the 8746B S-Parameter Test Set, and the 905A Coaxial Sliding Load for calibration runs. 

Analog-to-digital converters, together with a data logger and a T.I Terminal make it possible for 
the Network Analyzer to interface with a CRD main frame computer. Any given frequency band of 
the 8620A Sweep Generator may be swept with a step size of 0.1. 0.2, 0.5. 1. 2, or 5%. Steps of 1% 
have been chosen for this characterization work. This enables many data points to be taken so that 
spurious points are of no significance and confidence can be generated that just cannot be obtained 
when widely separated spot frequency measurements are taken. This "overview" approach is facilitated 
by a screen display of the calculated p and e values from which display a hard copy can be obtained. 
Smoothing by visual inspection and/or a drawing curve or straight edge completes the process. Alterna¬ 
tively, a tabulation of the measured values can be obtained. 

Several methods of measuring the complex permeability and permittivity are suitable for use with a 
network analyzer. Of these, the reflection technique (S n scattering coefficient) employing first a single¬ 
thickness. then a double-thickness sample placed in front of an RF short is the technique used here. 

Since the S-parameter terminals of the network analyzer have precision APC-7-type coaxial connec¬ 
tors, the sample holder was made to interface as smoothly as possible with this type connector. The 
precision outer (0.2756-inch, 7mm) and inner (0.1197-inch) conductors are available commercially and 
were used to construct the coaxial line sample holder. The dimensions of a toroidal sample are 0.275 
inch O.D., 0.120 inch I.D., and about 0.040 inch thick. (The sample thicknesses of the pair used for a 
particular set must be identical). An attempt is made to permit no more than 0.001 inch clearance 
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p = 0.4, /> EL = 55 X I0‘ 8 ftm 



Figure 11. Relative imaginary magnetic permeability, vs. frequency calculated for 
composites described in Figure 10. 


between the sample and its holder. In Figures 12-17 we show examples of these measurements for 
in the 6-12 GHz band for samples ONR #2, #4, #8. The first two are for powder with diame¬ 
ters in the range 20/am to 37 M m. loadings p = 0.4 and 0.2 respectively. The third of these samples is 
for powder with d < lO/am and p = 0.2. (cf. Table 3). The figures show the plotted data and the 
smooth curve used for subsequent evaluation. 

We are currently working to extend these results within the centimeter range. In addition we have 
requested our colleagues at GE-Re-Entry Systems Operation, Philadelphia, to carry out measurements 
at 35GHz. Initially it was planned in the Contract Proposal that all measurements above the 12 GHz or 
18 GHz range (specifically at 35 and 94 GHz) would be carried out at the Naval Research Laboratory. 
However, we were somewhat dismayed by their requirement of a 5 in. x 5 in. seamless square sample 
for the 35 GHz measurement in the arch method. By contrast, the equipment at GE-RSO uses a single 
1/4 in. diameter disc in a circular waveguide cavity. The measurements are in progress at the time of 
this writing. (This is another unanticipated reason for a subcontract which was not scheduled at the 
start of the program). 







0 


FREQUENCY GHz 


Figure 12. Relative real permeability measurements, fi t ‘ vs. frequency; range 6-12 GHz. 

Sample ONR #2, p = 0.413, sphere diameter 20 < d < 37j4m. Data curve for 
100 points and smoothed line. 



Figure 13. Relative imaginary permeability measurements, fi r " vs. f; same sample and fre¬ 
quency range as in Figure 12. 
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Figure 14. Relative real permeability measurements, n t ' vs. frequency; Range 6-12 GHz, 
Sample ONR #4, p = 0.200, sphere diameter 20 < d < 37/*m. Data and 
smoothed curve. 



IS. Relative imaginary permeability measurements, n” vs. f; same range and 
pie as in Figure 14. 
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Figure 16. Relative real permeability measurements, ft r ' vs. frequency; Range 6-12 GHz, 
Sample ONR #8, p = 0.208, sphere diameter d < 10um. 












11.2.4.1 Comparison of Permeability with Model Calculations 

From measurement curves of the type just presented we can read the smoothed curve values at a 
few (say 2 or 3) frequencies. A Gerber Variable Scale, Model TP007100B, is rather useful at this stage. 
These readings can be transposed onto graphs of the type shown in Figures 10, 11, i.e., the model cal¬ 
culations for vs. frequency for various volume loadings and various particle diameters. 

These calculations with the data for the seven samples measured in the range from 6 to 12 GHz are 
shown in Figure 18-25. These go in pairs, first p.,' then p.", each for the four p values (0.4 down to 
0.1). In two of the four p-value cases we have 2 or 3 particle size ranges. The agreement between 
measured data and calculated values is qualitatively quite good and even quantitatively rather acceptable 
in most examples. The experimental demonstration of the predicted and obviously non-trivial-in¬ 
magnitude behavior is an important step in the accommodation of artificial dielectric materials to mil¬ 
limeter wavelength applications. 

As we extend these measurements to higher (and lower) frequencies, the quality of the demonstra¬ 
tion will become more dramatic. At this point, the importance of this diversion in the research pro¬ 
gram has been clearly established. 

II.3 Microstructure of Artificial Dielectric Composites: Observations and Implications 

A significant task envisioned in the planning of this investigation, as shown in Table 1, and 
amplified in the Proposal Work Statement is to “Analyze the effects of microstructure control, e.g. 
porosity, passivation coating, packing, on the dielectric behavior of oxide-passivated metal-insulator 
composites.” This task was motivated by a basic question posed in Section I, i.e., the relationships 
between previously observed dielectric behavior and the simple Maxwell model, its range of validity out 
to higher [1] volume loadings than sustained by many other composites, and the observed deviations 
from the simple model behavior. 



Figure 18. Comparison of Experiment and Model Calculations, n,‘ vs. f, p = 0.4. 
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Figure 19. Comparison of Experiment and Model Calculations, /*/' vs. f; p = 0.4. 


We cannot complete a discussion of these questions and factors at this point, in part because the 
investigation is still in progress both experimentally and analytically, and in part because it depends on a 
more thorough exposition of theoretical treatments of composite dielectrics which comes later in this 
report. While this is a bit of a chicken and egg problem, we can make a start on the observations which 
are, or may be, possible and introduce their relevance in at least a sketchy way. 

II.3.1 The Insulating Coating 

In Section II. 1.5 we reported on the surface condition of the metal powder prepared by gas- 
atomization. From magnetic Curie temperature evidence, from Auger electron spectroscopy and from 
powder resistivity measurements we concluded that the atomized powder particles have an insulating 
coating. In several rough approximations we estimated the oxide coating thickness at about 0 1 to 
0.2/i.m. 

It would be rash to call this a passivating layer, however. This conclusion derives from the consid¬ 
erable reaction suffered by this alloy powder when composite preparation using an inorganic binder was 
attempted. That reaction was evidenced in the thermo-magnetic scan. 

The presence of the insulating coating is very significant for the electromagnetic behavior, especially 
for the dielectric properties. Loosely speaking, it prevents the individual metal particles from shorting 
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Figure 20. Comparison of Experiment and Model Calculations, p/ vs. f; p = 0.3. 


to each other. More elegantly, it means we have a correlated heterogeneous medium in which conduct¬ 
ing particles are completely surrounded by insulating material. This restricts the theoretical approach to 
describing the dielectric behavior in a favorable way. Volume loadings up to or approaching the packing 
limits can be achieved with coated particles (spheres) without the establishment of conducting paths, 
i.e.. percolation. 

With a coating on the conducting sphere, as well as a binder material, the resulting composite is a 
three-phase medium. However, it is possible to replace the coated sphere combination by a single, 
slightly modified, dielectric sphere. We shall return to this in a later section, but it suffices to say that, 
for the very thin coatings estimated herein, the coated metal sphere dielectric constant is sufficiently 
high compared to that of the binder that the properties of the resultant composite e' are unchanged. 

II.3.2 Calculation of Porosity 

An aspect of the microstructure of composites is the porosity, i.e., the volume fraction of holes or 
voids. It is straightforward to calculate this quantity using the densities of the separate components and 
their proportions ascertained from observable measurements. We refer back to definitions of symbols 
used in Section II. 1.6.3 for evaluation of volume loading. 

Consider a mass m c of the composite occupying volume v c . Then 

v c p c = m c = v c (p a f a v + p b f b v + p h f h v ), 
but, p h = 0, of course. Thus we have 

Pc - Pa f"a + Pb fb> or 


But with fh -> 1 -f a v - f b , we can write 


fb 


Pc - Pa (*a 
Pb 


24 










Figure 21. Comparison of Experiment and Model Calculations, /t r " vs. f; p = 0.3. 


- + p (— -1) 


Alternatively, using the formula for p in observables, we obtain 


The principal sources of error in this determination of porosity are the density values for the binder 
material and for the composite. By contrast, the alloy density and the various magnetizations (on a 
mass basis) are sufficiently well known or well measured to be much less important in estimating the 
possible errors in f h v . Therefore, we can write 


df h v 


a Pb 


a f h v 
a p c 


dp c 


(,_££, *!L _££.(!_£«) + £« ilc 


dp^ 

Pc 




Figure 22. Comparison of Experiment and Model Calculations, p/ vs. f; p = 0.2. 


When we substitute typical numbers into this formula we find: at the p = 0.4 level, a 2% error (uncer¬ 
tainty) in p b contributes about 0.01 to df£ and a 1% error in p c contributes a similar uncertainty (0.01) 
to dfh - , at the p = 0.1 level each of such a pair of uncertainties (0.01) in dfh follows from 1% errors in 
Pb or Pc- 

In Table 4 we present the result of the porosity calculation. Other characteristics of these samples 
such as volume loading and particle size were given in Table 3. The values of porosity are generally 
small. The scatter between positive and negative values (the latter are meaningless) say more about 
our uncertainties than about the composite samples. At this stage we conclude that the preparation 
method is sufficiently in control that porosity is not a significant factor for this set of samples to affect 
the electromagnetic parameters. 

II.3.3 Direct Observation 

Can we confirm the calculated conclusion that porosity in the polyurethane artificial dielectric com¬ 
posites is negligible? That would seem to be an easy task in these days when a scanning electron 
microscope (SEM) is often readily available. However, preparation of samples for observation is a spe¬ 
cial art which is needed in the case of these polymeric composites. 

Our Polymer Physics Unit was able to help out with observations on these composites with special 
equipment. In the normal course of events, an as-cut sample was lightly coated with Au and Pd (a few 
hundred Angstroms coat) to carry off electronic charge and heal. However, if the sample was micro- 
tomed at room temperature, the resulting SEM picture was a mess, i.e., a torn surface which told more 
about the cut than about the composite. Their new equipment (an Ultracut E) allowed them to carry 
out a cold microtome at about -150C using a glass knife. 

The SEM pictures shown in Figure 26 for one sample (ONR #2, p = 0.4, 20^m<d<37^m) are 
clean and interesting. There are faint but discernable circles which are cut(!) sections of the metal 
powder spheres. Their dimensions correspond to particle sizes used. Note that smaller sizes may be 
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Figure 23. Comparison of Experiment and Model Calculations, **," vs. f; p = 0.2. 


anticipated in the cut, and a few are distorted from the circular cross section. Their identity is 
confirmed by EDAX (Energy Dispersive Analysis by X-rays) which shows a high count for Ni and a 
lesser one for Cr on these circles in contrast to adjoining matrix. A number of holes are present. 
Some of these are obviously hemispherical “pull-outs.” Others are distorted but with excess lip material 
suggesting a pull-out. There are many regions densely filled with material of one kind or another. The 
latter observation is in accord with the calculations suggesting negligible porosity. At this volume load¬ 
ing the metal particle density is sufficiently high that one cannot get a feel for “irregular fluctuations" 
in the loading. 

It is a matter of concern for understanding the electromagnetic properties on a quantitative basis 
that one look for wide variations in local packing on some moderate length scale. Several workers have 
pointed out the potentially enormous effects that may arise from agglomeration with electric conduction 
contact [23] or wide variation in size distribution [24] intrinsically or as a consequence of 
agglomeration. In addition, as we discuss later herein, if a nominally dilute, low p-value composite has 
locally densely packed regions (of a size we hope to be able to investigate theoretically), compensated 
by regions nearly empty of metallic particles, without electric conduction paths in the densely packed re¬ 
gions, there can be some enhancement of electromagnetic property parameters. This arises from a 
strongly non-linear volume loading dependence of the property. 
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Figure 24. Comparison of Experiment and Model Calculations, n,' vs. f; p = 0.1. 

We have not yet made direct observations on composites of lower p value, but we plan to examine 
the p = 0.1 sample (ONR #5) in the near future. Also under consideration is the preparation of com¬ 
posite samples with still lower p values, for electromagnetic measurement and direct (SEM) observa¬ 
tion. 

II.4 Alternative Theoretical Approaches to Composite Dielectrics 

This section deals with an important task of the program. We have been quite successful in explor¬ 
ing a wide range of theoretical approaches from rather “classical” ones up to developments of the past 
half-dozen years, and in relating these to the materials physics of the applications of interest. This sub¬ 
ject attracted much attention in the nineteenth century (2) and yet it remains quite alive today as evi¬ 
denced by several recent topical conferences on the properties of heterogeneous media [18,25]. 

There is a basic distinction between several classes of models of such dielectric media. More 
modern theory follows the work of Bruggeman [2,26]. This is a symmetrical (two-component) effective 
medium model wherein particles of type 1 are mixed randomly in some proportion with those of type 2, 
so that each may be surrounded by particles of both types. If one type is conducting and the other 
insulating, there will be a percolation threshold composition at which conducting paths are established. 
The volume loading for this to occur is about 0.15 to 0.2. This model is often labeled EM A (eff. med. 
approx.). 

In another (older) class of model one assumes that particles of type 1 are completely surrounded by 
material of type 2. It is clearly necessary to specify whether metal has been added to a dielectric host or 
vice versa, i.e., the model is unsymmetrical. This class of model is associated with the names of 
Maxwell, Maxwell-Garnett, Clausius and Mossotti, and Lorentz and Lorenz. It too is really an effective 
medium (or molecular field) model. It is variously called by any of the names just cited, e.g., MGT, 
CM, etc. Real-life examples of this model are ordered (cubical) arrays of spheres in a different 
medium, or an otherwise correlated medium such as coated spheres in a random arrangement, with 
each metal particle surrounded by a concentric dielectric coating. In our prior work [1] we employed 
this model (under the name Maxwell model) without a critical evaluation and with moderate success. 
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Figure 25. Comparison of Experiment and Model Calculations, /* 


Table 4 

CALCULATED POROSITIES 
OF POLYURETHANE 
TEST COUPONS 


Sample # fg 










Figure 26. Scanning Electron Microscope views of cold-microtomed polymeric composite 
ONR #2 (p * 0.4, 20fim < d < 37^im). Top photo at 150 x; lower photo at 
500 x. 
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In recent years, rigorous calculations (and some measurements) of the ordered arrays have 
appeared [27,28,29,30] and they have relevance for the random coated array. Features of these theories 
are the significant effects of higher order multipoles beyond the dipole-dipole model and the (obvious, 
but sometimes ignored) maximum packing limits. Re-examination of the earlier data [1] confirm the 
applicability of these recent calculations, and deviations therefrom cast some light on the microstructure 
of the random coated arrays. 

II.4.1 The Maxwell (MGT.CM) Model 

An excellent derivation of the formulas for this model with emphasis on the concepts and approxi¬ 
mations is presented by Landauer [2], A contrasting experimentally oriented discussion on methods of 
approximating the behavior of two-phase (heterogeneous) systems has been developed by Mitoff [31], 

Mitoff’s report also calls attention to a very simple method described by Lichtenecker [32] for 
approximating the properties of two-phase systems, particularly when the properties of the two phases 
are not very different. Although the method is not fundamentally correct, it is easy to use and does 
agree with experimental results. 

Lichtenecker’s rule states that a plot of log k (k being the property value or “constant” for the 
composite) versus volume fraction p of phase 2 in phase 1 is a straight line between the end values. A 
graphical comparison with a correct treatment, such as Maxwell’s (discussed below) shows that the 
difference between the two is insignificant with respect to experimental accuracy where the constants of 
the two phases are within an order of magnitude of each other. 

Maxwell derived an equation which, at infinite dilution, exactly describes the k of a mixture of 
spheres with property constants ki and k 2 , (k 2 >k|). While he used k c and k d for the continuous and 
the discontinuous phases, Mitoff recasts this into two expressions: 

[ k| + 2k 2 - 2(1—p) (k 2 -k|) l 
k k2 [ ki + 2k 2 + (1-p) (k 2 —k!) | 
for phase 2 being the continuous phase; and 

f k 2 + 2k[ ~ 2p (k|—k 2 ) 1 

k k| [ k 2 + 2k, + p (k)—k 2 ) J 

for phase 1 continuous. It assumes that the discontinuous phase consists of spheres and that the 
spheres are separated enough not to disturb the lines of “flux” around a neighbor. It has been shown 
by experiment to be quite accurate up to 10% dilution and is a fair approximation in some cases even up 
to 50%. 

An alternative writing of the second of the above equations may be more familiar to some. Thus 
for the Maxwell, Maxwell-Garnett, Clausius-Mossotti model we also have 
3 p k| <k 2 —k,) 

k “ 1 + (1-p) k 2 + (2+p) k, 

When k 2 > >k, and p is small (low-concentration) this formula reduces toward 
k = k, (1 + 3 p). 

Although equations of the general type above have been available for many years, they have seen 
little use, probably because they are not the kind one keeps in one’s head for a quick approximation. 
Mitoff's report, however, shows that a very simple treatment applies when the two (separate) properties 
are very different. This is precisely the region for which the simple Lichtenecker method does not ap¬ 
ply. 


31 










For large differences, extending to a number of orders of magnitude, one adopts the procedure of 
plotting log k/kj versus volume fraction. From Mitoff, we obtain Figure 27. presenting the Maxwell 
curves for ki/ki ratios from 10 03 to 10 6 . The upper curves represent phase 2 continuous and the lower 
ones phase 1 continuous. The straight line, L, represents the Lichtenecker approximation. One notices 
that when kyk L is less than 10, the mathematically defensible Maxwell treatment is not much different 
from Lichtenecker’s. 

Mitoff goes a step further in observing that the limiting behaviors in Figure 27 for large differences 
(k 2 > > kj) are really straight lines with slopes that do not involve the property constant of the discon¬ 
tinuous phase. He obtains the relationships: 

A logio (k/k,) 


and 


A (1—p) 


A log l0 (k/k t ) 


-0.65 <k, = kJ 


These are shown graphically on Figure 28. So long as the k’s are very different, one expects the points 
to fall on these straight lines, especially at volume fractions close to the k c side. 


LOG k 2 /k, 


0.2 0.4 0.6 0.8 

VOLUME FRACTION, p 


1.0 


Figure 27. Maxwell (MGT, CM) Equations for Mixtures, on a Semilog Plot. Upper set is 
for Phase 2 continuous; lower set is for Phase 1 continuous; L is Lichtenecker's 
Approximation. (The number on each curve signifies the power of ten in the 
ratio k^k,.) (from Reference 31) 
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VOLUME LOADING, p 


Figure 28. MitofFs Graphical Approximations to the Maxwell (MGT, CM) Equations for 
k 2 > > k(. (from Reference 31) 

Mitoff compares this approximate method with the complete Maxwell equation and with 
Lichtenecker’s rule for a wide range of ratios of k^kp He shows that for one-order-of-magnitude 
difference or more the first or upper straight line for k 2 = k c is a fair approximation to Maxwell's result, 
and for two or more orders of magnitude difference the second or lower line for k| = k c is useful even 
out to rather high volume fractions of the discontinuous phase. At less than an order-of-magnitude 
difference, Lichtenecker's line is a good approximation to the Maxwell equation while the others are 
not. 

Inspection of Figure 27 shows that the Maxwell (MGT, CM) expressions continue analytically 
across the whole range of volume fraction, p. That the curves are shown dashed at extreme p values is 
a recognition of two points: one is the assumption of non-contact, translated as meaning a dilute sys¬ 
tem: second is the fact that the volume loading of uniformly sized spheres in random packing cannot 
exceed 0.63, and reaches this only under specially favorable conditions. We return to this point below. 
Nevertheless the Maxwell (MGT, CM) formula is useful, as is the Mitoff approximation when used 
with care. 

There are two applications of practical interest. The upper branch curves of Figure 27 (or their ana¬ 
lytic expression) successfully describe the effects of porosity in a host medium. The second application 
of special interest is that of a coaled sphere. Its relevance was noted in Section 11.3.1. An apparently 
direct application of the Maxwell (MGT. CM) expressions yields [33] the single dielectric constant 
2e, + e 2 + 2Q 3 (e 2 — e ( ) 

2e, +e 2 - Q J <e 2 - e,> 
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where 


Q = 1 - (t/a) 

with t the coating thickness and a the sphere radius. For e 2 »ei. medium 1 being the insulating coat, 
this reduces to 

1+2Q 3 
l—Q 3 

For our estimated coating thickness of 0.1/im, on a sphere of average diameter 23/xm, we find e cs /ei = 
115. For simplicity we take €j for the coating as that for Cr 2 0 3 , i.e., 11.1, from which e cs = 1276. This 
value for is still so much greater than that for typical binders that the model calculations are 
unaffected. 

II.4.2 The Bruggeman Effective Medium Approximation 

We described the initial elements of this EMA model in II.4 above. A visual image is presented in 
Figure 29 (after Landauer [2]). The two components are mixed randomly. The shaded crystallite of 
type 1 (of volume fraction, p) is surrounded by crystallites of both types. In the analysis [21, the sur¬ 
rounding (two-phase) material is imagined as a single medium of uniform properties. The homogene¬ 
ous dielectric constant e B in this case arises from the solution of a quadratic equation which, assuming 
spherical particles, is: 



P 


*B ~ < 1 
2« b + «i 


(1-p) 


g B ~ <=2 

2«B + «2 


0 


At quite low concentrations, p, the EMA and the Maxweil-MGT-CM give identical results. The EMA 
has the appealing aspect of treating particles of types 1 and 2 on an equivalent basis (hence a symmetri¬ 
cal theory). 



Figure 29, Schematic View for Bruggeman Effective Medium Model, (from Reference 2) 
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As we have discussed above, our artificial dielectric is one with a correlation between the two phases 
^ of the medium. This rules against the use of the EMA model and in favor of one for which the metal 

J is embedded in the insulating medium (thanks to the coating). 

An interesting and distinctive feature of the EMA is that is predicts a metal-insulator or percolation 
> transition. In this version it occurs at p = 1/3, while in other theories this value can be as low as p = 

0.15. (For quasi-two-dimensional mixtures, the values are somewhat higher). A good illustration of 
this prediction experimentally is shown in Figure 30 [34], For these samples the critical percolation frac- 

* tion (vol. loading) is about 0.20. This demonstration is a convincing argument against the EMA class 

H of model for the artificial dielectrics of this study. 

• II.4.3 Exact Calculations for Arrays of Spheres 

Concurrent with the interest in random heterogeneous media, is an attention to the equivalent prob¬ 
lem for ordered arrays of (conducting) elements. Lord Rayleigh in 1892 showed a significant advance 
_ beyond Maxwell’s work when he took into account induced octupole moments. The ordered array 

I lends itself to a more rigorous formulation of the problem. The subject received a technical boost with 

the suggestion of Kock [35] that artificial dielectrics (of the ordered type) could be used as microwave 
lenses. There has been a great deal of engineering interest, involving a variety of conductor 
shapes [36]. We have found particularly useful a paper by Kharadly and Jackson [27] both for the 
(lossless) effective permittivity calculations and measurements, and for discussions of the dielectric 
loss. Most recently, with modern computers, one can carry out numerical calculations of these prob- 
< lems to arbitrary precision limits which were totally unattainable heretofore [28-30,37,38], 
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As noted above, the higher order multipoles (beyond the dipole terms of the Maxwell/CM theories) 
have an important effect as the spherical conductors get closer at increased p. A measure of the success 
of the calculations is that the calculated permittivities exhibit an anomaly at the critical volume loading 
where the spheres just touch, in contrast with the Maxwell (MGT, CM) models. The conceptual drama 
of this feature shows in the simplicity of a confirmatory experiment for the body centered cubic lat¬ 
tice [30], It suffices to construct only a one-eighth portion of a unit cell using conducting and insulat¬ 
ing plates plus two one-eighth sections of metal ball bearings! The metal-insulator "anomaly” or 
transition when they touch is thus lattice dependent. This follows from the 3-dimensional packing lim¬ 
its for hard, uniform spheres, as follows: 

simple cubic p = it 1 6 = 0.52 

body centered cubic p = n V3/8 = 0.68 

face centered cubic p = ir >1216 = 0.74 
or hexagonal closed packed 

In Figure 31 we reproduce the results of the calculations, the upper from Doyle [28,37] and the 
lower from the Australian group [29,30,38], The Maxwell-CM curve in the upper portion agrees with 
the exact calculations out to p^ 0.35 for SC and to p =0.5 for BCC and FCC. This provides the 
justification for the reasonable success in our prior work. It was the right type of model and a not-bad 
approximation. On the other hand, the Maxwell/CM model fails to diverge until p tends toward unity. 
The line marked R is the Rayleigh calculation for the SC lattice (corrected [28]) showing some mul¬ 
tipole effects but failing to diverge at the contact limit. In addition experimental data are included for 
the SC lattice which are in line with the exact results. 

The calculations of the Australian group in the lower portion of Figure 31 agree in detail with those 
of Doyle, since the methods of attack were equivalent. They have added another curve estimated by 
scaling arguments and intended to apply to a random or disordered array. The empirical packing limit in 
this case (as noted above) is p = 0.63, although that limit is tricky to reach and values nearer to p = 
0.60 are more likely. They even express this curve with a closed approximate formula as follows: 

e = 0.05398—2.06751n(0.63—p). 

II.5 Dielectric Measurements and Analysis 

In discussing the magnetic dipole effects ( induced! in artificial dielectrics (Section II.2 above) it was 
important to examine both the real and imaginary parts of the permeability. Following the infrared 
studies of small metallic particles [17,19] it was noted that electric dipole effects contribute significantly 
to absorption in that frequency range only when the metal particle size is below about 50A. The size 
dependence of the magnetic dipole loss effect is one factor here. The other is that for appropriate 
models of the dielectric behavior of a metal [17,19,27] the frequency region for significant electric loss 
is governed by the intrinsic electron plasma frequency of the metal and the electron relaxation time. 
These correspond to frequencies in the ultraviolet and visible ranges, respectively, clearly far away from 
the microwave and millimeter frequency ranges of present interest. 

Therefore, to the extent that we can consider our composite binders to approximate perfect (loss¬ 
less) dielectrics, we need focus only on the real part, e’, in our measurements. From an operational 
point of view, however, we do not ignore the measurement results which include e" values. It is 
important to examine the measurement output for c" in that “unsmall values” are often a sign of 
erratic operation of the measurement process. 

In what follows we shall revisit the dielectric behavior of composites included in our prior work with 
a view to reconciling this with the more modern theories. Finally we present and discuss the e’ data 
obtained on the composites prepared for this study. 








VOLUME RATIO, p 



Figure 31. Dielectric Constants for Cubic Lattices as calculated in Ref. (28) (upper) and 
Refs. (29, 30) (lower). 


II.5.1 Reanalysis of Prior c' Results 

In our prior study [1] we obtained dielectric measurement data on a composite system of magnetic 
metal particles in the same polyurethane binder as employed herein. The particles had an insulating 
coating (c.f. Section II.1.5.1 and Table 2). The composites studied had measured volume loading rang¬ 
ing from about p = 0.25 to p = 0.50. It is instructive to re-examine these data in the light of 
Section II.4, inasmuch as the prior analysis had considered only the Maxwell (MGT-CM) model. 

In Figure 32 we display that «’ data set along with a number of comparative curves. The various 
symbols used for the data set refer to series with differing particle sizes and size ranges. Despite some 
experimental scatter, the family of data describe a general curve. The solid straight line passing through 
the lower p-value data is the MitofT approximation to the Maxwell model. (Recall from Figure 28 that 
its slope is pre-determined.) In the limit of p = 0. this line comes close to, but deviates a bit from the 
«' data point measured for the binder alone. We more-or-less ignored this deviation at the time of that 
work and we do not now have an explanation for it. 










VOL FRACTION, p 


Figure 32. Measured Dielectric Constant Data from Prior Study; also Comparative Model 
Curves for Cubic Lattices (SC, BCC, FCC) and an Estimated Random Array 
Curve (R). 

The upward deviation of the data from the Mitoff approximation line is a more prominent feature. 
In the prior study we tried to attribute the deviation to voids, without any analysis. That attempted ex¬ 
planation contradicts the effects of voids as we presently understand them. However, we do not have 
to look far for a qualitative explanation. The offset family of curves in Figure 32 are obtained from 
Figure 31, the calculated curves of «' vs. p for cubic lattices, including higher order multipole effects. 
(The dotted curve is the Mitoff line, again). The general feature of the upward curvature of the e' vs. 
p data is clearly a consequence of the multipole interations and the limits to packing. 

From a quantitative point of view the explanation leaves room for comment. The data curve rises 
slightly more rapidly than the calculated curve for the simple cubic (SC) lattice. It is markedly away 
from the curve estimated for a random distribution. In discussions with Prof. W.T. Doyle of Dart¬ 
mouth College (References 28, 37), he suggested an attractive explanation. 

The suggestion is based on the marked non-linearity of the e' vs. p curve. In contrast to the lattice 
arrays of particles, one could have fluctuations in volume loading in a randomly packed composite. 
Indeed fluctuations are the norm in random systems but we imagine that on some spatial scale it can 
make a difference. If, for example, a nominal p = 0.50 sample were comprised of equal parts of p = 
0.4 and p = 0.6 loadings, the ratio f/«i would be 5.20 instead of the estimated random ratio of 4.27. 
(The separate values of «/«i for p = 0.4, 0.6 are taken from the estimated curve/formula in 
Section II.4.3). 
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We remarked on this aspect of the random composites in Section II.3.3 in connection with direct ob¬ 
servation of the composite in a SEM. For the one p = 0.4 sample shown in Figure 26, such effects are 
not easily detected. By examining some more dilute samples we hope to find evidence for such (non¬ 
contacting) clumping. 

In addition there may be some interesting questions about the spatial scale of clumping needed to 
produce the electromagnetic effects hypothesized. We plan to pursue these and related points in a 
collaboration/consultation with Professor Doyle. 

II.5.2 Measurements on Non-Magnetic Composites 

The artificial dielectric composite samples described in Section II. 1.6 were measured by the tech¬ 
niques of Section II.2.4. Major emphasis at this stage has been on the frequency range 6-12 GHz. 
Extensions to lower and higher frequencies are underway, including measurement at 35 GHz. Future 
plans call for measurements at 94 GHz. (Permeability results on these samples were reported in 
Section II.2.4.1). In Figure 33, we show examples of these measurements for 3 samples, p = 0.4, 0.3, 
0.2 of the Ni 92 Cr 8 powder in polyurethane with 20/im < d < 37^m. (Samples ONR #2, #3, #4). 

An overview of the e' vs. p data is presented in Figure 34, along with a few model curves. Of the 
latter, it is useful to see how closely the Mitoff approximation line comes to matching the exact simple 
cubic array dielectric curve. In considering the measured data, picked at 10.5 GHz where the noise is 
minimal, it is clear that we have a problem. The data could be fit reasonably with a Mitoff line, but 
only if we ignore the binder-only point (taken anew for this study). 



FREQUENCY, GHz 


Figure 33. Relative Real Permittivity, «/, vs. Frequency; Range 6-12 GHz. Samples ONR 
#2, 3, 4. 
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Figure 34. Overview of Relative Permittivity, c' vs. p Various particle sizes; N^Crg in 
Polyurethane Composites. 

It seems rather far to stretch to attribute this upward shift to clumping in the way just described 
(See 11.5.1). At this stage, with limited data, there does not seem to be any significant variation with 
particle size. We would not expect any under ideal conditions. 

We shall continue this investigation in depth in order to resolve the problem. We plan to make (by 
the same techniques) samples with lower p-values. The independent measurements currently underway 
at 35GHz may confirm or deny the accuracy of the techniques used. It is additionally puzzling that the 
prior data shown in Figure 32 exhibited this effect so modestly that we tended to ignore it. 

This concludes our report on the group of Tasks numbered 1, 2 and 3 in Table 1. 

III. TASK 5: EXPLORE ALTERNATIVE PASSIVATION COATINGS 

The system of interest in the prior work [1] has been called an oxide- passivated metal-insulator 
dielectric. The oxide coating was obtained by a classified technique. In the present study, scheduled 
for the second year, there is a task of exploring alternative or improved coatings. We have made a start 
on this task and details are presented in a classified appendix to this Annual Report. 





IV. TASK 6: QUANTIFY SKIN-DEPTH EFFECTS IN PASSIVATED 
ARTIFICIAL DIELECTRICS 

From our paraphrasing of the Statement of Work, this task can be posed under the question - Has 
the artificial (magneto-) dielectric system been optimized in the centimeter wavelength regime? 
Achieving such a goal carries with it the possibility of predicting the particle sizes needed for millimeter 
wavelength applications. 

It was our intent to conduct this portion largely empirically, and then deduce the guiding rules which 
would facilitate prediction. There is however, some reason to feel that a reasonably analytic approach 
may be possible, in parallel with the experiments. 

The main challenge in this task is in the area of fine metal powder preparation, and processing. 

IV.1 Spark Erosion Preparation of Fine Powders 

We have made a promising start toward the preparation of powders for high-frequency optimization 
of the artificial magneto-dielectric system. We are exploiting the potential of the spark erosion tech¬ 
nique to prepare suitable crystalline alloy particles. This technique has enjoyed renewed attention of 
late with important contributions from our colleagues. 

A schematic of the spark erosion apparatus is shown in Figure 35. The process of spark ero¬ 
sion [39] results from maintaining a repetitive pulsed spark discharge between electrodes immersed in a 
dielectric fluid. The spark discharge produces highly localized superheated regions in the electrodes. 
Small molten droplets and vaporized materials are ejected from these heated regions into the dielectric 
fluid, where they are rapidly cooled into particles. This process was used to prepare colloids at least 
50 years ago. In more recent times it is the basis of electric discharge machining. The choice of dielec¬ 
tric fluids can range from organics through aqueous and cryogenic liquids. 



Figure 35. Spark Erosion Apparatus. The electrodes and the coarse charge are of the same 
composition. The larger particles (>0.5>tm) which form from erosion of the 
charge and electrodes settle through the screen. 
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Previous workers have noted that the particle yield could have a bimodal distribution of diameters. 
Smaller particles with average diameters about 50A tend to remain in colloidal dispersion. “Larger” 
particles, with diameters between 0.5 and 25/*m, settle out. Our colleagues have recently examined the 
smaller ones with a view toward their application as metallic ferrofluid particles (instead of the usual 
magnetite, Fe 3 0 4 , ferrofluid material) [39]. 

In the present work our immediate goal is to prepare and separate powders, and then evaluate com¬ 
posites of magnetic alloy particles in the less than lpt range. This will be a continuation of prior work 
which studied artificial dielectric material of various particle sizes down to a range of 1 to 4/*m diameter 
in an insulating binder [1], Secondarily this task addresses the question of just how effective the spark 
erosion method might be as an alternative for general fine powder alloy production. 

A large number of experiments have been undertaken using electrodes and charge of several suit¬ 
able alloys. The dielectric fluids have included liquid argon, liquid nitrogen, ethyl-alcohol, water and, 
most recently, mixtures of the last two. Another significant variable in these tests is the duration of the 
spark pulse, e.g., short at 3/xs and long at 50/as. 

With the aqueous fluids, there is a distinct tendency for the particle size distribution to shift toward 
smaller particles witn shorter spark discharge pulses. There have been major (>50%) yields in the 
d< 10/Ltm range. This is extremely interesting when compared with existing commercial techniques. 

The technological potential inherent in cheap dielectric fluids, such as water or water-based mix¬ 
tures, is also very attractive. It permits a freer rein in design conception when one considers possible 
scale-up. 

On the other hand there is a greater tendency for oxidative processes to accompany the particle for¬ 
mation in an aqueous fluid. We are just beginning to sort out these effects and hence a discussion of 
detailed results will be deferred to a future report. 

IV. 2 Possible Analytic Approach to Skin-Depth Effects in Artificial Dielectrics 

The goals of optimizing the system (i.e., the metal particle size) in the centimeter wavelength range 
and forecasting the properties at higher frequencies were posed originally as semi-empirical tasks. The 
progress reported in Section II.2. on understanding and demonstrating induced magnetic dipole and par¬ 
ticle size effects has changed our view of the possibilities. The development at that point was for a 
strictly non-magnetic artificial dielectric. The present task carries the implied goal of involving the 
magneto-dielectric system of our prior study [1], 

In that earlier work we examined, in a very simple model, some eddy current screening effects on 
the intrinsic magnetic permeability (/ju"). We did not however, incorporate this into a global model. 
In the course of the present work we encountered an old study which presents an overall attack on fre¬ 
quency dependence of such heterogeneous systems. On the other hand, it omits the ferromagnetic 
resonance considerations which formed an important part of our prior study. We feel that there is an 
interesting potential here in revisiting that older work and marrying it to the additional concepts just 
noted. We take this as a supplementary goal in the coming period of investigation. 

V. CONCLUSIONS AND FUTURE PLANS 

Inasmuch as we are just one-third of the way into this extended study, the matter of “conclusions” 
is only partly that, and partly an accounting of achievements en route to goals. It is thus appropriate to 
follow such a status review with an enumeration of plans for the coming period. 

V.l Conclusions 

During the first phase of the study we have developed an interesting artificial dielectric. It is 
comprised of spherical alloy powder particles of the composition Ni 92 Cr g prepared by gas atomization 
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with liquid substrate quenching which produces an insulating oxide coating on each particle. The com¬ 
posites fabricated with the metal powder and polyurethane are ferromagnetic only below room tempera¬ 
ture, a feature which allows for easy determination of volume loading. Before mixing with the binder, 
the particles are classified by size at <10/*m, +10— 20/tm, and +20—37^m. 

The decision to pay close attention to powder size followed from the realization that examination of 
induced magnetic dipole (eddy current) effects would be as important as the dielectric studies originally 
planned. Calculations of complex magnetic permeability, adapted from the literature, showed that the 
dominant parameter is the ratio of particle radius to skin depth. The sizes chosen and frequencies of 
study (centimeter to millimeter wavelength ranges) cover the most interesting behavior anticipated. 
Measurements so far, mainly in the 6-12 GHz band, agree with calculations within about 20%. The 
results have applications consequences for understanding the properties of metallic magneto-dielectrics 
under extreme conditions of temperature or frequency, as well as for their frequency dependence under 
“normal” operating conditions. 

Exploration of theoretical developments in the dielectric behavior of heterogeneous media has suc¬ 
cessfully addressed one of our opening questions, i.e., why our prior (and present) material systems fol¬ 
low so well the simple Maxwell (Maxwell-Garnett, Clausius and Mossotti) theoretical model. The in¬ 
herently isolated metallic particles represent a correlated system even when loaded randomly in a 
binder, in contrast to the percolation-prone Bruggeman model exemplified in many experimental metal- 
insulator composites. Modern theories which incorporate higher order electric multipole effects and 
geometrical packing limits have been validated previously only (to our knowledge) by ordered lattice¬ 
like arrays of spherical particles. Our random composites also exhibit such properties, albeit with devia¬ 
tions probably reflecting on their microstructure. 

We started work on alternative passivation coatings, of interest scientifically and technologically. We 
have also made good progress toward the preparation of very fine powders by the spark erosion method, 
with the goal of an empirical optimization of the metallic magneto-dielectric system in the centimeter 
wavelength range. The initial results on particle size distribution are quite interesting with future tech¬ 
nological potential. The use of simple dielectric fluids including aqueous mixtures or water alone also 
holds technological promise. 

V.2 Future Plans 

We are quite interested in exploiting possible relationships between our combined experimental and 
theoretical success with the complex permeability of the non-magnetic artificial composites and the 
recent decade of results and calculations on anomalous absorption in the far infrared by fine metal par¬ 
ticle composites. The anomalies of the latter studies represent a disagreement between theory and 
experiment of several orders of magnitude, in contrast to our reasonable success. To be sure, our fre¬ 
quencies and sizes differ from those of the infra-red studies by about a factor of ten (each) but that is 
not so far away considering the contrasting behaviors. We hope to investigate the connections, and 
perhaps to explore the possibility of collaboration with one of the academic groups working in the infra¬ 
red. Essential to this is some extension of our own work into the millimeter wavelength range, a task 
which is a bit behind schedule, in part because of changes in scope. 

The dielectric property phase of the study holds several interesting challenges. One is to explore 
experimentally the possible clustering in randomly packed composites toward understanding the devia¬ 
tions in the e' vs. p behavior. Some theoretical approaches to dielectric properties as related to the 
dimensions of fluctuations in packing may be undertaken. We hope to work with Prof. W T Doyle of 
Dartmouth on this aspect. Perhaps related to these aspects is the problem of apparently poor correla¬ 
tion in our new artificial dielectric system between the trend of «’ vs. p and the p = 0 limit (binder 
only) measured value. We need to understand this discrepancy. 

Another problem area is the fabrication of a composite system using a second binder of noticeably 
different 

We intend to pursue the alternative coatings investigations as outlined. 
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There is some excitement in the future as regards the goal of optimizing the magneto-dielectric sys¬ 
tem in the centimeter wavelength range and predicting the properties in the millimeter range. Experi¬ 
mentally this excitement centers on the research and technology potentials associated with fine powder 
production by the spark erosion method. Analytically we are cautiously hopeful that our progress in 
understanding the induced magnetic dipole effects can be extended and broadened to be a quantitative 
predictive tool for metallic magneto-dielectric systems. We plan to push this as well. 
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